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We correct one error which appeared in the paper Asteroseismological constraints on the pulsating planetary 
nebula nucleus (PG1159-type) RX J 2111. 1+3412, A&A, 461, 1095 (2007), by A. H. Corsico, L. G. Althaus, M. 
M. Miller Bertolami & K. Werner (DOI: 10.1051/0004-6361:20066452), regarding the age of expansion of the 
planetary nebula of RX J2117. 1+3412 (£pn)- The correct value is ipN ~ 2.71 x 10 4 yr, instead of the published 
value (tpN ~ 5.43 x 10 4 yr). Note that the evolutionary timescale of the central star of RX J2117. 1+3412 is 
t* ~ 2.5 x 10 4 yr. Thus, by considering the correct expansion age, we definitively solve the discrepancy between 
the RX J2117. 1+3412 evolutionary timescale and the size of the nebula. The notable agreement between £pn 
and i* reinforces the validity of our asterosesimological model for RX J2117. 1+3412. 



1 



Astronomy <fc Astrophysics manuscript no. text6452 


© ESO 2008 


February 5, 2008 





Asteroseismological constraints on the pulsating planetary nebula 
nucleus (PG1159-type) RX J2117.1+3412 

A. H. Corsico 1 ' 2 *, L. G. Althaus 1 ' 2 *, M. M. Miller Bertolami 1 ' 2 ' 3 ** and K. Werner 4 

1 Facultad de Ciencias Astronomicas y Geofisicas, Universidad Nacional de La Plata, Paseo del Bosque S/N, (1900) La 
Plata, Argentina. 

C^- ' 2 Instituto de Astrofi'sica La Plata, IALP, CONICET-UNLP 
<^^> 3 Max-Planck-Institut fur Astrophysik, Garching, Germany 

4 Institut fur Astronomie und Astrophysik, Universitat Tubingen, Sand 1, 72076 Tubingen, Germany 
O^l ■ e-mail: acorsico , althaus ,mmiller@f caglp .unlp. edu. ar ; werner@astro.uni-tuebingen.de 

5_l ' 

■ Received; accepted 



(N 
> 
O 
(N 
^- 
O 

o 



- 1—^ 

X 



ABSTRACT 



Aims. We present asteroseismological inferences on RXJ2117. 1+3412, the hottest known pulsating PG1159 star. Our 
results are based on full PG1159 evolutionary models recently presented by Miller Bertolami & Althaus (2006). 
Methods. Wc performed extensive computations of adiabatic g-mode pulsation periods on PG1159 evolutionary models 
with stellar masses ranging from 0.530 to O.741M0. PG1159 stellar models are extracted from the complete evolution of 
progenitor stars started from the ZAMS, through the thermally pulsing AGB and born-again phases to the domain of 
the PG1159 stars. We constrained the stellar mass of RX J2117. 1+3412 by comparing the observed period spacing with 
the asymptotic period spacing and with the average of the computed period spacings. We also employed the individual 
observed periods to find a representative seismological model for RX J2117. 1+3412. 

Results. We derive a stellar mass M, ~ 0.56 — O.57M0 from the period spacing data alone. In addition, we found 
a best-fit model representative for RX J2117. 1+3412 with an effective temperature T c g = 163 400 K, a stellar mass 
A/„ = O.565M0, and a surface gravity log g = 6.61. The derived stellar luminosity and radius are log(L*/L0) = 3.36 and 
\og(R* / Rq) — —1.23, respectively, and the He-rich envelope thickness is M em = O.O2M0. We derive a seismic distance 
d ~ 452 pc and a linear size of the planetary nebula Dpn ~ 1.72 pc. These inferences seem to solve the discrepancy 
between the RXJ2117. 1+3412 evolutionary timescale and the size of the nebula. All of the seismological tools we use 
concur to the conclusion that RX J2117. 1+3412 must have a stellar mass M* ~ O.565M0, much in agreement with 
recent asteroseismology studies and in clear conflict with the predictions of spectroscopy plus evolutionary tracks. 
C/5 ■ 

Key words, stars: evolution — stars: interiors — stars: oscillations — stars: variables: other (GW Virginis) — white 
dwarfs 



$H 1 1. Introduction MacDonald 2003, Althaus et al. 2005, and Miller Bertolami 
■ ■ ^ • , et al. 2006) or a late helium thermal pulse (LTP) that takes 
Pulsating PG1159 stars (or GW Vir stars) are very hot place durfng the post _ AGB evolution when hydrogen burn- 
hydrogen-deficient post-Asymptotic Giant Branch (AGB) ing ig stm actiye (gee B15cker 20 01 for references). During 
stars with surface layers rich m helium, carbon and oxy- a yLTp episode; most of the hydrogen-rich envelope of the 
gen (Werner & Herwig 2006). They exhibit multiperiodic star is burnt in the helium nash -driven convection zone, 
luminosity variations with periods m the range 300 - 3000 while in a LTP the hydrogen deficiency is the result of 
seconds, attributable to nonradial pulsation g-modes. Some a dilution episode. In both cases, the star returns rapidly 
GW Vir stars are still embedded m a planetary neb- back tQ the AQB and finally into the domain of high effec _ 
ula (commonly called PNNVs; Planetary Nebula Nucleus tiye tem perature as a hydrogen-deficient, quiescent helium- 
Variable), and are characterized by much higher lummos- burning object . This helium burning phase continues until 
ity than the "naked" GW Vir stars (those without neb- the hot whi te-dwarf domain is reached, 
ulae). PG1159 stars are thought to be the evolutionary 

link between post-AGB stars and most of the hydrogen- RXJ2117. 1+3412 is the hottest known pulsating 
deficient white dwarfs. These stars are believed to be the PG1159 star. It was discovered as an X-ray source in the 
result of a born again episode triggered either by a very ROSAT sky survey (Motch et al. 1993). The star is em- 
late helium thermal pulse (VLTP) occurring in a hot white bedded in an old, extended and diffuse planetary neb- 
dwarf shortly after hydrogen burning has almost ceased (see ula, which at the distance of 1400 pc (Appleton et al. 
Schonberner 1979 and more recently Herwig 2001, Lawlor & 1993) makes it the largest presently known with a diam- 
eter of 5.3 pc. This implies an expansion age much in ex- 
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mospheres, T c g = 170 ± 10 kK and log g = 6+Q2, and a 
surface composition of He/C/0= 47.5/23.8/6.2 (by num- 
ber of atoms). Furthermore, a strong mass- loss rate of 
M « 10 -7 — 10 _8 M Q /yr has been measured (Koesterke 
et al. 1998; Koesterke & Werner 1998). The central star of 
RX J2117. 1+3412 was discovered to be a pulsator indepen- 
dently by Watson (1992) and Vauclair et al. (1993), with 
pulsation periods that are intermediate between those of 
PNNVs and naked GW Vir stars. The effective temperature 
of the central star defines the blue edge of the observational 
PG1159 instability strip. Vauclair et al. (2002) (hereinafter 
VEA02) published an amazing richness of asteroseismolog- 
ical results on RXJ2117. 1+3412 from a multisite photo- 
metric campaign that included two Whole Earth Telescope 
(WET, Nather et al. 1990) runs. Without actually perform- 
ing pulsation computations, they were still able to derive 
fundamental quantities for RX J2117. 1+3412, such as the 
stellar mass, the helium-rich envelope mass fraction, rota- 
tion rate, luminosity, and distance. 

It is important to note that at the time the analysis 
of VEA02 was performed, no realistic PG1159 evolution- 
ary models suitable for interpreting high-luminosity, high- 
temperature transition objects such as RX J2117. 1+3412 
were available. In fact, while the 0.7M Q evolutionary se- 
quence of Wood & Faulkner (1986) available at that 
time fits the location of RX J2117. 1+3412 in the logg - 
logT c ff plane, its predicted surface abundance composi- 
tion is not representative of the observed abundances in 
RX J2117.1+3412. More importantly, the Wood & Faulkner 
(1986) models do not constitute appropriate structures for 
PG1159 models that have been evolved through a born- 
again episode. On the other hand, the evolutionary tracks of 
Gautschy (1997), with surface abundances more appropri- 
ate to PG1159 stars, do not have luminosities high enough 
to be compatible with RXJ2117. 1+3412. In view of the 
lack of realistic stellar models for a full asteroseismologi- 
cal analysis, VEA02 were forced to use the scaling relation 
between period spacing and mass of Kawaler & Bradley 
(1994) to infer the stellar mass of RX J2117.1+3412. It is 
worth noting that this relation was derived for GW Vir 
stars characterized by much lower luminosities than those 
of the PNNVs. 

One of the results of VEA02 is that the asteroseismologi- 
cal mass they derive for RXJ2117. 1+3412 is about O.56M , 
in serious disagreement with the mass derived from spec- 
troscopical analysis (0.72 M & ). This fact has motivated us 
to undertake the present investigation taking full advan- 
tage of the new generation of PG1159 evolutionary models 
recently developed by Miller Bertolami & Althaus (2006). 
In fact, these authors have computed realistic PG1159 evo- 
lutionary sequences for different stellar masses by consider- 
ing the complete evolution of progenitor stars. They have 
followed in detail all of the evolutionary phases prior to 
the formation of PG1159 stars, particularly the born-again 
stage. We believe that the evolutionary models presented 
by Miller Bertolami & Althaus (2006) represent a solid the- 
oretical background to analyze the evolutionary and pulsa- 
tional status of hot PG1159 stars like RX J2117. 1+3412, a 
transition object for which stellar models calculated with 
artificial evolutionary procedures should be taken with cau- 
tion. 

In this work we perform an asteroseismological analy- 
sis of RXJ2117. 1+3412 based on the PG1159 evolutionary 
models of Miller Bertolami & Althaus (2006). Emphasis is 



placed on deriving pertinent structural parameters for this 
star by using its observed period spectrum. In particular, 
we derive the stellar mass for this object from three different 
seismological approaches. Additionally, the implications of 
our results for the evolutionary status of RX J2117. 1+3412 
and for its planetary nebula are discussed. The paper is or- 
ganized as follows: in the next Section we briefly describe 
the input physics of the PG1159 evolutionary models and 
the pulsational code we employ. In <j3]we derive the stellar 
mass of RXJ2117. 1+3412 by means of the observed period 
spacing. In Section|3]we derive structural parameters of this 
star by employing the individual observed periods. In Sj5]we 
discuss the discrepancy we found between the seismological 
and the spectroscopic mass of RXJ2117. 1+3412, and in fJB] 
we discuss the thickness of the outer envelope of the star as 
predicted by our seismological model. The asteroseismolog- 
ical distance to RXJ2117. 1+3412 and the implications of 
our analysis for its planetary nebula are provided in Section 
[7] Finally, in Sect. [8] we summarize our main results and 
make some concluding remarks. 

2. Evolutionary models and numerical tools 

The pulsation analysis presented in this work relies on stel- 
lar models that take into account the complete evolution 
of PG1159 progenitor stars. The pulsation stability proper- 
ties of these PG1159 models have recently been analyzed by 
Corsico et al. (2006). The stages for the formation and evo- 
lution of PG1159 stars were computed with the LPCODE 
evolutionary code, which is described in detail in Althaus 
et al. (2005). Briefly, LPCODE uses OPAL radiative opac- 
ities (including carbon- and oxygen-rich mixtures) from 
the compilation of Iglesias & Rogers (1996), complemented 
at the low-temperature regime with the molecular opaci- 
ties of Alexander & Ferguson (1994) (with solar metallic- 
ity). Chemical changes are performed via a time-dependent 
scheme that simultaneously treats nuclear evolution and 
mixing processes due to convection, salt finger and over- 
shooting. Convective overshooting is treated as an expo- 
nentially decaying diffusive process above and below any 
convective region. 

Specifically, the background of stellar models were ex- 
tracted from the PG1159 evolutionary calculations recently 
presented in Miller Bertolami & Althaus (2006) (0.53, 
0.542, 0.565, 0.609 and 0.664 M sequences), Althaus et al. 
(2005) (0.589 M Q ) and Corsico et al. (2006) (0.741 M ). 
The range of initial masses of these models on the ZAMS 
is 1 — 3.75-Mq. We refer the reader to those works for de- 
tails. All of the sequences were followed from the ZAMS 
(assuming solar metallicities) through the thermally puls- 
ing and mass-loss phases on the AGB. After experiencing 
several thermal pulses, the progenitors departed from the 
AGB and evolved towards high effective temperatures. The 
mass-loss rate during the AGB was arbitrarily kept fixed, as 
to obtain a final helium shell flash during the early white 
dwarf cooling phase (VLTP). As a result, evolution pro- 
ceeds through the born-again stage, which brings the rem- 
nants back to the giant domain. After this episode, most 
of the hydrogen content of the stars was burnt, and the 
hydrogen-deficient remnants evolved at constant luminos- 
ity to the domain of PG1159 stars with surface chemical 
composition rich in helium, carbon and oxygen. The evo- 
lutionary tracks in the logT e ff — logg plane are displayed 
in Fig. [TJ For the 1-Mq sequence two different AGB evolu- 
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Fig. 1. The Miller Bertolami & Althaus (2006) PG1159 
full evolutionary tracks in the logT c ff — log g plane, 
corresponding to stellar masses of (from right to left): 
0.530,0.542,0.565,0.589,0.609, 0.664, and 0.741M Q . Thick 
(thin) curves correspond to stages before (after) the models 
reach their highest effective temperature. The black circle 
with the error box is the location of RXJ2117. 1+3412 ac- 
cording to spectroscopy (T c ff = 170 ± 10 kK and \ogg — 
®-o 2) an d the black square is the location of the star as pre- 
dicted by our asteroseismological analysis (Jl]). The blue 
(hot) and red (cool) boundaries of the theoretical dipole 
(£ = 1) instability domain according to Corsico et al. (2006) 
are depicted with thin dashed lines. 

tions were considered, with different mass-loss rates as to 
obtain different numbers of thermal pulses and, eventually, 
two different remnant masses of 0.530 and 0.542M Q . Mass- 
loss episodes after the VLTP were not considered in the 
PG1159 evolutionary sequences we employed here, despite 
the fact that mass-loss episodes for PG1159 stars have been 
reported (see Werner & Herwig 2006). 

We computed adiabatic pulsation periods with the help 
of the numerical code we employed in our previous works 
(see Corsico & Althaus 2006 for details). We analyzed 
about 3000 PG1159 models covering a wide range of ef- 
fective temperatures (5.4 ^ log(T c g) ^ 4.8) and luminosi- 
ties (0 <; log(L»/L©) <; 4.2), and a range of stellar masses 
(0.530 < M*/M Q < 0.741). For each model we computed 
£ = 1 g-mode periodfQ in the range 50 <; II fc <; 3600 s, 
which comfortably covers the period spectrum observed in 
RXJ2117. 1+3412 (690 - 1200 s). 

3. Mass determination from the observed period 
spacing 

Here, we constrain the stellar mass of RXJ2117. 1+3412 by 
comparing the asymptotic period spacing (AIT|) and the 

1 Quadrupole (I = 2) modes have not been detected in 
RXJ2117.1+3412 (VEA02). 



average of the computed period spacings (AII/s) with the 
observed period spacing (An j3. These methods take full 
advantage of the fact that the period spacing of PG1159 
pulsators depends primarily on the stellar mass (Kawaler 
& Bradley 1994; Corsico & Althaus 2006). 

The asymptotic period spacing and the average of the 
computed period spacings for £ = 1 modes as a function 
of the effective temperature are displayed in Figs. [2] and 
[31 respectively, for different stellar masses. Also shown in 
these diagrams is the location of RX J2117. 1+3412, with 
T off = 170 ± 10 kK (Werner et al. 1996) and AlT 3 = 
21.618 ± 0.008 s (VEA02). Here, AITf = Tl /^£{£+ 1), 
where n = 2tt 2 [J^ (N/^dr}' 1 , being N the Brunt- Vaisala 

frequency (Tassoul et al. 1990). The quantity Allfc, on the 
other hand, is assessed by averaging the computed forward 
period spacings (AHk = Uk+i — life) in the range of the 
observed periods in RXJ2117. 1+3412. 

From a comparison between Aid and Ald| ("first 

method" ), we obtain a stellar mass of M* = O.5681q'oo3-^0- 
The quoted uncertainties in the value of M* come from 
the errors in the spectroscopic determination of the effec- 
tive temperature^- Note that, since the M* = O.565M0 
evolutionary track does not reach the spectroscopic T e g of 
RXJ2117.1+3412 (see Figs. Q] and [2]), we are forced to ex- 
trapolate the value of A1T| for that effective temperature 
(170 kK) and also for its upper limit (180 kK) to derive the 
stellar mass. 

In the same way, we get M* = O.56OIq;q*3-M0 if we com- 
pare All and All^ ( "second method" ) . Again, we need to 
extrapolate the values of Allfc for the spectroscopic value 
of T c g and its upper limit in order to derive the stellar mass 
and its lower limit. 

It is evident that both estimations of the stellar mass 
nearly agree Q. The slightly higher value of M* as derived 
from Allf is due to that usually the asymptotic period spac- 
ing is a bit larger than the average of the computed period 
spacings (see Corsico & Althaus 2006). The first method to 
derive the stellar mass is somewhat less realistic than the 
second one, because the asymptotic predictions are, in prin- 
ciple, only valid for chemically homogeneous stellar models, 
while our PG1159 models are indeed chemically stratified 
(see g6]). 

VEA02 derive a stellar mass of M* = O.56i^M by 
extrapolating IIo from the formula of Kawaler & Bradley 
(1994). It is interesting to note that, despite that extrapola- 
tion being based on a relation which is valid for luminosities 
between 1.6 < log(I/*/L Q ) < 2.8 (much lower than the lu- 
minosity of RX J2117.1+3412), the mass derived by VEA02 

2 Note that the vast majority of asteroseismic studies are 
based on the asymptotic period spacing to infer the stellar mass 
of GW Vir pulsators. 

3 Our uncertainties are different from those quoted in VEA02 
which are due entirely to the fact that they lack evolutionary 
calculations appropriate for RX J21f 7.1+3412. 

4 Note that in both derivations of M» we have assumed that 
RX J21f 7.1+3412 is evolving towards higher effective tempera- 
tures, before it reaches the evolutionary knee in the logT c ff — 
logg diagram. In other words, we have used the thick portions 
of the curves in Figs. [5]and|3]to infer the stellar mass. Assuming 
instead that the star is on the beginning of its white dwarf cool- 
ing phase, then the thin portions of the curves should have been 
considered. In that case we would have obtained values of the 
stellar mass only slightly different from those quoted in the text. 
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Fig. 2. The dipole (£ = 1) asymptotic period spacing 
(All!) f° r different stellar masses in terms of the effec- 
tive temperature. Numbers at the left of each curve de- 
note the stellar masses (in solar units). Thick (thin) curves 
correspond to stages before (after) the models reach their 
highest effective temperature. The small black circle with 
the error bar depicts the location of RX J2117. 1+3412 
(T eff = 170 ± 10 kK and AIT 3 = 21.618 ± 0.008 s). 

is in complete agreement with our predictions, although our 
uncertainties are much smaller. 

Finally, we note that our inferred stellar mass value of 
M* « 0.56 — O.57M , is in serious conflict with the value 
M* = 0.72+0 ^2-^0 as derived from spectroscopy coupled 
to evolutionary tracka^l (Miller Bertolami & Althaus 2006). 
This will be discussed in more detail below. 

4. Constraints from the individual observed periods 

In this approach we seek a pulsation model that 
best matches the individual pulsation periods of 
RX J2117. 1+3412. The goodness of the match between the 
theoretical 1=1 pulsation periods (life) and the observed 
periods (IIP) is measured by means of a quality function 
defined as X 2 (M*,T eff ) = £™ =1 m i n [(np - U k f]/n, where 
n (= 20) is the number of observed periods (first column 
in Table [J). The PG 1159 model that shows the lowest 
value of x 2 Wln be adopted as the "best-fit" model. 

We evaluate the function x 2 (Af*, T e g) for stellar masses 
0.530,0.542,0.565,0.589,0.609,0.664, and 0.741M o , 
whereas for the effective temperature we employed a much 
more finer grid (AT off = 10 - 30 K). The function X 2 
in terms of the effective temperature for different stellar 
masses is shown in the mosaic of Fig. [4] We found only 
two important minima of comparable magnitude, one of 
them for the model with A/* = 0.565M Q and T cff « 163.4 

5 The quoted uncertainties in the spectroscopic value of M* 
have their origin in the uncertainties in the determination of 
logg and T e a, e.g., we neglect possible uncertainties in the evo- 
lutionary computations. 
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Fig. 3. Same as FigJ2l but for the average of the computed 
period spacings (Allfe). 



kK (panel c) and the other one for the model with 
M* = 0.589A/© and T eS w 105 kK (panel d). Another 
minimum, albeit less pronounced, is also encountered for 
the model with M* = 0.565Af Q at T cS « 123 kK. We adopt 
the first model as the best-fit model because it has an 
effective temperature relatively close to that suggested by 
spectroscopy. A detailed comparison of the observed m = 
periods of RX J2117. 1+3412 with the theoretical periods 
of the best-fit model is given in Table [TJ The high quality 
of our period fit is quantitatively reflected by the average 
of the absolute period differences, <5IT = \^-i\)/ n i 
where <5IT = 11^ — life , and by the root-mean-square resid- 
ual, a sn . = v^(E |<5H,| 2 )/n. We obtain 5n 4 = 1.08 s and 
a s ^ = 1.34 s. The quality of our fit for RXJ2117. 1+3412 
is substantially better than those achieved by Kawaler & 
Bradley (1994) and Corsico & Althaus (2006) (STh = 1.19 
s and Wi = 1.79 s, respectively) for PG 1159-035. 

Table Q] also shows the linear growth rates (j]k) of the 
fitted pulsation modes (fifth column), computed with the 
nonadiabatic pulsation code described in Corsico et al. 
(2006). A positive value of T)k implies pulsational instabil- 
ity. It is interesting to note that the domain of unstable- 
mode periods (757 <; life <; 1167 s) of our best-fit model 
nearly coincides with the range of the observed periods for 
RX J2117. 1+3412, although our calculations predict the ob- 
served modes with k = 30, 31, 32 and 53 to be pulsationally 
stables, and conversely, modes such as k = 38,41,48 and 
52, which are unstable according to our predictions, are not 
observed at all in RXJ2117. 1+3412. 

The last column in Table Q] shows the theoretical rate 
of period change of the fitted pulsation modes. Our cal- 
culations predict all pulsation periods to decrease with 
time (life < 0), an effect attributable to the rapid con- 
traction experienced by the star in its evolution to the 
blue (see Fig. [T]) . Unfortunately, an unambiguous II mea- 
surement for any mode of RX J2117. 1+3412 is currently 
very difficult because of the large amplitude variations of 
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T eff = 170 000 K, log g = 6.0 (spectroscopic result) 
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Fig. 5. Detail of a HST spectrum of RX J2117. 1+3412 comprising a Civ line and the location of the T e g sensitive 
Ov 1371 A line. The Ov line is not detected, placing a lower limit to T g. Upper panel: Comparison with a synthetic 
spectrum from a model with T c g=170 kK and log g =6.0, the best-fit parameters from optical and UV spectroscopy. 
At T o ff<170 kK the model predicts a detectable Ov line. Lower panel: Comparison with four synthetic spectra with 
T e ff=160, 170, 180, 190 kK. The O v line gets weaker with increasing T c g . The surface gravity is fixed at log g =6.6, the 
asteroseismic result. At this higher gravity, T e g must exceed 180 kK in order to have non-detectable O V line, significantly 
above the asteroseismic result of T c ff=163 kK. 



the modes (VEA02). On the basis of considerations of 
the mass-loss rate and the helium-rich envelope mass in 
RXJ2117.1+3412, VEA02 give a rough estimation of the 
rate of period change of |fl| ~ 3 — 24 x 10~ n s/s, much in 
agreement with our predictions. 



5. Discrepancy between spectroscopic and seismic 
mass 

The main features of our best-fit model are summa- 
rized in Table [2E where we also include the param- 
eters of RXJ2117. 1+3412 from other published stud- 
ies. Note that the effective temperature of the best-fit 
model is compatible at la with the spectroscopically in- 
ferred value. According to our results, the temperature of 
RXJ2117. 1+3412 is somewhat lower than previous estima- 
tions, but it is still the hottest known PG1159 pulsator. 
On the other hand, the best-fit model has a stellar mass 
of M* = 0.565M Q , very close to the values found from 
the period spacing data (see JJ]). Thus, the linear pulsa- 
tion theory strongly suggests that RX J2117. 1+3412 is an 
average-mass PG1159 star with M* ~ 0.56 — 0.57M Q . This 
is clearly at odds with the spectroscopic-evolutionary pre- 
dictions, that indicate a rather high mass of 0.70 — 0.72M Q 
for this star (Werner & Herwig 2006; Miller Bertolami 
& Althaus 2006). A mass discrepancy is also encountered 



Errors in T e g and \og(L* / L@) are estimated from the width 
of the minimum in the function \ 2 vs T e s and \og(L, / Lq), re- 
spectively; the error in the stellar mass comes from the grid res- 
olution in M». Errors in the other quantities are derived from 
these values. 



with other PG1159 pulsators (PG 1159-035, PG 1707+427, 
PG 2131+066, PG 0122+200). In those cases the seismic 
masses are always higher than the spectroscopic masses (see 
Corsico et al. 2006), however, the differences are a factor of 
three smaller (0.04 - O.O6M ) than for RX J2117. 1+3412. 

The mass discrepancy for RX J2117. 1+3412 is reflected 
by the difference between the surface gravity of our pulsa- 
tion model and the value predicted by spectroscopy. In fact, 
for our best-fit model we have log<? = 6.61, large in excess 
as compared with the derived value of logy = 6.0^2 by 
using NLTE model atmospheres. Hence, the location of the 
star in the log T e g — log g plane shifts, according to our seis- 
mological predictions, downwards well beyond the limits of 
the observational error box (see Fig. [1} . 

It is obvious that the seismic results on T e s & n d 
log g are difficult to reconcile with results of model atmo- 
sphere analyses. This is demonstrated with Fig.O which 
shows a HST spectrum of RXJ2117. 1+3412 in the vicinity 
of the O v 1371 A line. This feature is very temperature 
sensitive and is therefore used as a strategic line to fix T e ff 
in PG1159 stars. It is not detectable in RXJ2117. 1+3412 
which puts a tight lower limit to T e ff. The lower panel 
of Fig. [5] shows a series of model profiles with increasing 
T e ff, whereas the gravity is kept fixed at the asteroseismic 
value log g = 6.6. The model atmospheres were calculated 
with the Tubingen Model Atmosphere Package (see, for de- 
tails, Werner et al. 2003). The model with T eff =160 kK (a 
value close to the asteroseismic result of 163 kK) produces 
a strong O v line, clearly contradicting the observation. A 
model temperature of 190 kK is necessary to make the line 
undetectably weak. In the upper panel we show the model 
with the spectroscopically determined best-fit parameters 
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minimum seen in panel c (arrow). The vertical dashed line 
is the spectroscopic T off of RXJ2117. 1+3412 (170 kK) and 
the grey zone depicts its uncertainties (±10 kK). 



(T e ff =170 kK, log g =6.0). In the analysis of Werner et al. 
(1996) we regarded this T e ff as a lower limit but excluded 
temperatures higher than 180 kK because all He n and C iv 
lines in the optical and UV become too shallow. This dis- 
crepancy between spectroscopic and asteroseismic results 
calls for a comprehensive re-analysis taking into account all 
available spectroscopic material collected since the Werner 
et al. (1996) analysis. 

VEA02 infer the stellar luminosity and radius of 
RX J2117. 1+3412 on the basis of their mass estimation and 
the effective temperature and gravity measured by spec- 



Table 1. Observed m = periods (n°) of 
RXJ2117.1+3412 (taken from VEA02), theoretical 
I = 1 periods (life), period differences (<5IT), radial orders 
(k), normalized growth rates (%), and rates of period 
change (life) of the best-fit model. 
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troscopy. Thus, not surprisingly, their results point to a 



= larger luminosity and radius (logL*/i0 = 4.05 



lGg(JZ./iZo) = - 
model (logT»/T Q 



91 + 010 

u - yi -0.15 



-0.23 
-0.32 



and 



as compared with our best-fit 
3.36 and log(i?*/i? Q ) = -1.23). 



6. Helium envelope thickness 

Another important parameter is the thickness of the outer 
envelope (M env ) of RXJ2117. 1+3412. We define M cnv as 
the mass above the chemical discontinuity between the He- 
rich envelope and the C/O core (see Fig. [5]). Our best-fit 
model has M cnv = O^MgQ. VEA02 claim that M cnv is 
0.0073 — 0.044Mq, which is in line with our results. It is 
important to note that VEA02 derive these values of M onv 
on the basis of the results published by Kawaler & Bradley 
(1994) for PG 1159-035, extrapolated to the range of pa- 
rameters of RX J2117. 1+3412. Specifically, they use a rela- 
tion connecting the fractional mass of the helium-rich en- 
velope with the effective temperature and the trapping cy- 
cle for PG 1159-035 and RXJ2117. 1+3412. The underlying 
assumption is that the mode-trapping features of PG1159 
stars are inflicted mainly by the outer chemical interface of 
He/C/O, and thus, the minima seen in the period-spacing 
distribution are associated with modes trapped in the en- 
velopes. However, realistic evolutionary calculations predict 



7 We define the "effective" location of the He/C/O chemical 
interface as the position of the layer characterized by an helium 
abundance of Xhc ~ 0.19. In our best-fit model this corresponds 
at Mr ~ 0.544M© (see the vertical line in Fig. [6}, where the 
contribution of the He/C/O chemical interface to the Brunt- 
Vaisala, frequency is largest. 
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mode-trapping properties of PG1159 models to be primar- 
ily fixed by the stepped structure of the C/O chemical pro- 
file at the core left by prior extra mixing episodes (Fig. 
EJ, the outer He/C/O interface playing a secondary role 
(see Corsico & Althaus 2005, 2006). For instance, the ob- 
served periods at 799.50 s (k = 35) and 1058.03 s (k = 47), 
which correspond to modes trapped in the envelope of 
RX J2117.1+3412 according VEA02, are modes "confined" 
to the core region in our best-fit model and thus not related 
to the thickness of the helium-rich outer envelope. Thus, the 
procedure from which VEA02 derived the thickness of the 
helium-rich outer envelope of RXJ2117. 1+3412 should be 
taken with some caution. 

The employment of full PG1159 evolutionary models 
allows us to place constraints on the amount of stellar 
mass eroded by winds. We can do this by adopting typical 
mass loss rates for the different stages after the VLTP. For 
the short superwind phase, that lasts ~ 800 yr and takes 
place at log T c g < 3.8, we adopt 10~ 5 M Q /yr as observed 
at Sakurai's Object (Tyne et al. 2002). For the evolution 
before the star reaches T ff~ 100 000 K (this stage takes 
^2000 yr) we adopt typical [WR] planetary nebula nuclei 
mass loss rates of 10 -6 M Q /yr (Koesterke 2001). Finally, 
for the PG1159 stage we can adopt the observed value of 
4 x 10~ 8 M Q /yr (during the 12000 yr until the sequence 
reaches the best fit model). Under these assumptions we 
estimate the possible mass lost by the star after the VLTP 
to be 10~ 2 Mq, which is half of the mass of the He-rich 
envelope. Thus the remaining envelope will be still thick 
enough for it not to affect the results of the present work. 
Even more, as the evolution of the star would be acceler- 



ated by these mass loss episodes, the previous estimation 
represents an upper limit under the adopted mass loss rates. 

7. The seismic distance to RX J2117. 1+3412 and 
the implications for its planetary nebula 

In closing, we infer the seismic distance of 
RXJ2117.1+3412 from the Earth. First, we consider 
the flux predicted by a NLTE model atmosphere with 
T ff = 160 kK and logg = 6.6 integrated through the 
spectral response of the V filter. We obtain a bolomet- 
ric correction BC= —7.954 and an absolute magnitude 
M v = 4.43. Because the proximity of RXJ2117. 1+3412 
to the galactic plane, we must account for the interstel- 
lar absorption, Ay. In line with VEA02, we follow the 
interstellar extinction model developed by Chen et al. 
(1998). We compute the seismic distance d according to 
the well-known relation: 



logd 



1 



Ay(d)], 



(1) 



where the apparent magnitude is m v = 13.16 ± 0.01 
(Motch et al. 1993). The interstellar absorption Ay(d) 
varies non linearly with the distance and also depends on 
the Galactic latitude (b) and longitude (£). For the equa- 
torial coordinates of RXJ2117. 1+3412 (Epoch B2000.00, 
a = 21 h 17 m 7 S .60, S = +34° 12' 22". 0) the correspond- 
ing Galactic coordinates are b = —10° 24' 32". 04 and 
£ = 80° 21' 10". 8. We solve Eq. © iteratively and ob- 
tain a distance d — 452 pc and an interstellar extinction 
Ay = 0.45, substantially lower than the estimations of 
VEA02 (d = 760lJ3° pc and Ay = 0.86). This is due to 
the fact that the luminosity of our model (2.3 x 10 3 L Q ) is 
considerably lower than the luminosity adopted by those 
authors (1.1 x 10 4 L o ). 

Our calculations predict a parallax of ~ 2.21 mas and 
a linear diameter of the planetary nebula of ~ 1.72 pc 
(5.3 x 10 13 km). If we assume a mean expansion velocity for 
a typical planetary nebula with a WR-type central star of 
approximately 31 km/s (Gorny & Stasinska 1995), then the 
expansion age for the nebula would be ipN ~ 5.43 x 10 4 yr. 
On the other hand, our computations predict an evolution- 
ary age for the best-fit VLTP model of about t* ps 2.5 x 10 4 
yr and possible larger if the VLTP would have occurred 
later on the hot white dwarf cooling branch than in our 
sequencfl This age is well consistent (to a factor less than 
two) with the expansion age of the nebula. Thus, our re- 
sults seem to solve the apparent inconsistency between the 
evolutionary timescale of RX J2117. 1+3412 and the size 
of the nebula (see Appleton et al. 1993 for details). The 
agreement between t* and ipN according to the present cal- 
culations reinforces the validity of our asteroseismological 
model and may suggest that RX J2117. 1+3412 have under- 
gone a VLTP episode, a fact that would be at variance with 
the lack of 14 N in its atmosphere. 



8. Summary and conclusions 

In this paper we carried out an asteroseismological anal- 
ysis of the hot pulsating PG1159 star RX J2117. 1+3412, 



We mention that our O.565M0 sequence predicts a time in- 
terval between the occurrence of the VLTP and the present state 
of RXJ2117.1+3412 of w 1.5 x 10 4 yr. 
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Table 2. The main characteristics of the central star and the planetary nebula of RX J2117. 1+3412. The second column 
corresponds to spectroscopic results, whereas the third and fourth columns present results from the pulsation study of 
VEA02 and from the asteroseismological model of this work, respectively. 



Quantity 


Spectroscopy 


Pulsations 
(VEA02) 


Asteroseismology 
(This work) 


An [s] 




21.618 + 0.008 


21.707 


Tea [kK] 


170 ± 10 (a) 




163.4+3^ 


M, [Ms] 
logs [cm/s 2 ] 
log(WL ) 
\og{R t /R Q ) 

M env [Mq] 


„ 79 +0.15(b) 

a n +0.3(a) 
O.U_ Q 2 

3.95 ± 0.5 (c) 


0.56l° ;S 
4.051°;! 

-0 Q1+ 010 
u - 91 -0.15 

0.0073 - 0.044 


565+ 024 

U.OOO_ 23 
a ci +0.11 

6-61-0.07 
3.36 ± 0.04 

i oq+0. 046 
1 - ZO -0.025 

0.02 ± 0.006 




39 55 Of! 

\J.OcJ) XJ.OOj u.uu 




39 39 22 


BC [mag] 






-7.954t° ;°j 


My [mag] 
A v [mag] 
d [pc] 
7r [mas] 

t. w 


140 q+700(c) 

71 + a4 ° 

u - ' i -0.24 

2 - 5 x 10 4(d) 


0.86 

760 + 23 ° 
' DU -235 
i OO+0.59 
J - oz -0.30 

> 1.3 x 10 5 


44 n + 0.12 

n 4t -+0.036 
u - 4o -0.012 

452 +4 ? 

o 91+O.12 
z - zl -0.20 

> 2.5 x 10 4 


D PN [pc] 
ipN [yr] 


r n + 2.7(d) 
d -°-1.9 

> 1.5 x 10 5(d) 


2.9 + 0.9 

~ 9.1 x 10 4 (*) 


1 7 9+0.18 
L - ' -0.09 

5.43t°;g x 10 4 



References: (a) Werner et al. (1996); (b) Miller Bertolami & Althaus (2006); (c) Motch et al. (1993); (d) Appleton et al. (1993). 
(*) This value is inferred from the linear size of the planetary nebula derived by VEA02 and a mean expansion velocity of the 

nebula of 31 km/s. 



a 5-mode pulsator with properties intermediate between 
young planetary nebula central stars and PG1159 stars 
without planetary nebula. Our analysis is based on the full 
PG1159 evolutionary models of Miller Bertolami & Althaus 
(2006). These models represent a solid basis to analyze the 
evolutionary and pulsational status of hot PG1159 stars 
like RXJ2117. 1+3412, a transition object for which stellar 
models extracted from artificial evolutionary procedures are 
not appropriate. We first made good use of the fact that the 
period spacing of variable PG1159 stars is mainly a function 
of the stellar mass, and derived a value of M* ~ 0.568M Q 
by comparing Aid with the asymptotic period spacing of 
our models. We also compared Aid with the computed 
period spacing averaged over the period range observed in 
RXJ2117. 1+3412, and derived a value of M* - O.56OM . 
Note that in both derivations of the stellar mass we made 
use of the spectroscopic constraint that the effective tem- 
perature of the star should be of ~ 170 kK. 

Next, we adopted a less conservative approach in which 
the individual pulsation periods alone naturally lead to 
an "asteroseismological" PG1159 model that is assumed 
to be representative of RXJ2117. 1+3412. The period fit 
was made on a grid of PG1159 models with a quite fine 
resolution in effective temperature (AT e g ~ 10 — 30 K) al- 
though admittedly coarse in stellar mass (AM* ~ 0.012 — 
O.O77M ). The match between the computed dipole pulsa- 
tion periods of the best-fit model and the observed periods 
in RX J2117. 1+3412 is excellent, with an average of the ab- 
solute period differences of 1.08 s and a root-mean-square 
residual of 1.34 s. It is worth noting also that the domain of 
unstable-mode periods of the best-fit model nearly matches 
the range of the observed periods in RXJ2117. 1+3412, al- 
though there is no obvious correlation between the mag- 
nitude of the linear growth rates and the observed mode 
amplitudes. 



Interestingly enough, the mass of the best-fit model 
(M* = 0.565Mq) agrees with the mass derived from the ob- 
served period spacing, but it is ~ 25% lower than the value 
M* = 0.72M Q derived from spectroscopy coupled with our 
evolutionary tracks. 

Other characteristics of the best-fit model are summa- 
rized in Table [5J In particular, the effective temperature is 
compatible at la with the spectroscopic value. At variance 
with this, the surface gravity of the best-fit model is sub- 
stantially higher than that given by spectroscopy. We also 
infer the "seismic distance" of RX J2117. 1+3412 and ob- 
tain a distance d ~ 452 pc, which places the star markedly 
closer to the Earth than thought hitherto (d ~ 760 — 1400 
pc, Vauclair et al. 2002 and Motch et al. 1993). The derived 
distance implies a linear size of the nebula of 1.72 pc which 
implies an expansion age of w 5.43 x 10 4 yr. This age is 
substantially lower than assumed hitherto (Appleton et al. 
1993) and in better agreement with the times predicted by 
evolutionary models. This result reinforces the correctness 
of our asteroseismological model for RXJ2117. 1+3412. 

Finally, our computations predict a temporal period 
drift for RX J2117.1+3412 between -3 x lO" 11 s/s and 
— 1.4 x 10~ 10 s/s. The negative values of Id reflect the fact 
that our best-fit model is still rapidly contracting on its 
evolutionary road towards higher effective temperatures in 
the HR diagram. Unfortunately, the amplitude variability 
of the observed modes in RX J2117. 1+3412 precludes any 
measurement of II for the moment, thus hindering any test 
of our prediction. 

The results of the asteroseismological analysis car- 
ried out in this work strongly suggest that the mass of 
RXJ2117. 1+3412 is considerably lower than suggested by 
spectroscopy coupled to evolutionary tracks (Werner & 
Herwig 2006; Miller Bertolami & Althaus 2006). This seri- 
ous disagreement is also seen in other pulsating PG1159 
stars like PG 1159-035, PG 1707+427, PG 2131+066 and 
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PG 0122+200, although in those cases the seismic masses Wood, P. R., & Faulkner, D. J. 1986, ApJ, 307, 659 
are always higher than the spectroscopic masses (Corsico 
et al. 2006). The discrepancy could be attributed to a num- 
ber of factors. On the observational side, possible errors in 
the spectroscopic determinations of g and T e ff, in particu- 
lar, as this study suggests in the case of RXJ2117. 1+3412, 
in the determination of g. In this respect, detailed astero- 
seismological fits to other pulsating PG1159 stars would 
be valuable. On the theoretical front, different PG1159 
evolutionary tracks could result from different dredge-up 
and/or mass loss history in the AGB progenitor evolution 
(Werner & Herwig 2006). But in preliminary simulations 
we have found that neither third dredge up efficiency nor 
TP-AGB lifetimes play a determining role in the location 
of PG1159 tracks. It remains to be seen if other assump- 
tions in the microphysics such as radiative and conductive 
opacities and/or equation of state may be playing a role in 
the location of post- AGB tracks. 
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